Phonon Heat Conduction in Corrugated Silicon Nanowires Below the Casimir Limit 
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The thermal conductance of straight and corrugated monocrystalline silicon nanowires has been 
measured between 0.3 K and 5 K. It is demonstrated that the corrugation strongly reduces the 
thermal transport by reducing the mean free path of the phonons. The experimental averaged 
mean free path is remarkably smaller than the smaller diameter of the nanowire, evidencing a 
phonon thermal transport reduced below the Casimir limit. Monte Carlo simulations highlight that 
this effect can be attributed to significant multiple scattering of ballistic phonons occuring on the 
corrugated surfaces. This result suggests an original approach to transforming a monocrystalline 
material into a phonon glass. 



PACS numbers: 63.22.-m, 63.20.dd, 66.70.Df, 65.40.-b 

In the past decade, the thermal and thermodynamic 
properties of nanostructured materials or low dimen- 
sional materials has attracted growing interest QUI], es- 
pecially for applications in thermoelectrics [5H7] . Recent 
experiments and theoretical research have shown the po- 
tential for engineering phonons at the nanoscale, for in- 
stance, the reduction in phonon thermal conductivity in 
nanostructured thin films, nanowires [2| I5HTT]. Overall, 
size effects on the heat conduction have been widely stud- 
ied, demonstrating effects on the phonon mean free path 
(MFP) [T2J[13], on the dispersion relation [14], or on the 
transmission coefficient [8j[T5]. In general, reducing the 
phonon thermal transport below the Casimir limit, where 
the phonon mean free path is limited by boundary scat- 
tering, remains an experimental challenge [T6HT9] . All 
these recent experiments have involved the measurement 
of systems with nanoscale sizes, which has motivated the 
development of specific fine thermal experiments [20H24] . 

In reduced size systems, surface scattering of phonons 
plays a crucial role [25] [26]. As a result, effects linked 
to surface roughness have been at the core of numer- 
ous studies showi ng a significant impact on the phonon 
transport especially at low temperatures [TOJ [TTJ I27U32] . 
Recent experiments have demonstrated very low thermal 
conductances of highly rough silicon nanowires [7] fT8] . 
far below the amorphous limit. Hence, the possibility 
of creating a phonon glass type material from a perfect 
monocrystal seems possible by playing on the geometry 
of the surface. These nanoengineered surfaces may intro- 
duce phonon backscattering or multiple scattering trap- 
ping phonons in a lattice of cavities [33]. In this respect, 
phononic crystals (a regular modulation of the geometry 
of a phonon waveguide) have attracted significant atten- 
tion in the past decade as a promising means to reduce 
phonon transport [34H38] . Despite numerous theoretical 
works, no measurement has been yet reported demon- 
strating phonon transport below the Casimir limit in pe- 



riodic (corrugated) phonon conductor. 

Here, we report experimental evidence of a strong re- 
duction in the phonon mean free path (MFP), below 
the Casimir limit, in silicon corrugated nanowires. The 
phonon flux is estimated through the measurement of 
the thermal conductance of straight and corrugated sus- 
pended nanowires. The modulated surfaces induce a se- 
vere reduction in the MFP, which becomes significantly 
smaller than the minimum diameter. Experimental data 
for the mean free path are compared with predictions 
based on the Monte Carlo technique, which involves mul- 
tiple scattering. 

The thermal conductance of suspended monocrys- 
talline silicon nanowires has been measured for tempera- 
tures ranging between 0.3 K and 5 K. The Si nanowires 
(straight and corrugated) are made by high resolution 
e-beam lithography from intrinsic silicon on insulator. 
The minimal cross section of the corrugated nanowire is 
150 mux 200 nm and the maximal one is 250 mux 200 nm 
[39] . The semi- length of the nanowires is either 4.4 /xm 
or 5 fim [40] and the roughness is estimated to be of a 
few nanometers [10] . Scanning Electron Microscope im- 
ages of the measured nanowires are shown in Fig. [I] The 
periodicity of the corrugation is 200 nm. This length as 
well as the width is still larger than the dominant phonon 
wavelength, estimated to be on the order of 65 nm at 1 K, 
and 6.5 nm at 10K : Xdom = 4 25 fc^ r > wnere ^ * s the 
Planck constant, v s refers to the speed of sound, which 
was set to 6500 m/s (the mean value of the speed of sound 
of the transverse and longitudinal acoustic phonons) , and 
fc# is the Boltzmann constant [41] [42]. between 1 K and 
10 K, all the sample dimensions are at least two times 
bigger than Xdom excluding significant effect of phonon 
confinement on the thermal properties. 

The experiments were performed by using the 3a; 
method [43 -45 . This measurement technique has al- 
ready been explained in detail [TO] [45] . A highly sensitive 
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FIG. 1: SEM images of the straight a) and the corrugated b) 
nanowires; c) corresponds to the top view of the corrugated 
nanowire. The scale bars correspond to a) 2 /mi, b) 2 /im and 
c) 300 nm. 



thermometer in niobium nitride is deposited on top of the 
intrinsic suspended silicon nanowires [46]. An RMS cur- 
rent smaller than 1 nA is used to create a temperature 
oscillation in the thermometer resistance (on the order 
of 600 kOhm at 1 K). This leads to a dissipated power 
ranging from 10 -15 Watt to 10 -13 Watt; a temperature 
oscillation amplitude of 3 to 10 mK is thus generated in 
the temperature range of 0.3 K to 6 K. 

The measured thermal conductance values are spread 
over three orders of magnitude between 10 -13 W/K to 
10 -10 W/K. All the measurements were made at a fre- 
quency below 5 Hz, where the nanowires are in the qua- 
sistatic thermal regime [45]. In this frequency range, the 
thermal signal is only sensitive to the thermal conduc- 
tance; the experiment has a resolution estimated to be 
a few percentage points with a high energy resolution 
(down to the atto Joule (10 -18 J)) [47] . The parasitic 
thermal path through the thermometer can be estimated 
by the Wiedemann-Franz law giving a thermal conduc- 
tance of 4.10 -14 W/K at IK, much below the thermal 
conductance of the nanowires. 

In Fig. [2] we show the thermal conductance data ver- 
sus temperature for the corrugated (samples CI to C4) 
and the straight nanowire. Table [T] gathers all experi- 
mental data and nanowire characteristics. The values 
obtained for the straight nanowire serve as a reference 
when analysing the thermal properties of the corrugated 
nanowires. In order to compare the nanowire thermal 
conductances to each other, they will be normalized with 
respect to the one predicted by the Casimir model. In the 
Casimir theory adapted by Ziman and coworkers [25] [48] , 



FIG. 2: (colour online) Thermal conductance versus temper- 
ature for a straight nanowire and four corrugated nanowires 
in the log-log scale. 



the thermal conductance is given by: 



K Ca 



3.2xlO d 



2ir 2 k% 



(2/3) 



e X wAcas _ 3 



(i) 

^Cas — l-12\/e x w is the phonon Casimir MFP for 
a rectangular shape phonon conductor, e refers to the 
thickness and w to the average width of the nanowire, L 
being its length [9]. This formula is obtained when the 
mean free path of the phonons is limited by the cross 
section of the wire, i.e., by boundary scattering. 

As the temperature decreases, the dominant phonon 
wavelength actually grows larger than the roughness 
height and specular reflections become more probable. 
To take this trend into account, Ziman proposed replac- 
ing the Casimir mean free path by an effective MFP A e ff 

EiHH]: 
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where p is a parameter describing the probability for a 
phonon to be specularly reflected when hitting the sur- 
faces. The value of this parameter will depend on both 
the wavelength of the phonons and the roughness of the 
wires. In this model, the phonon boundary scattering 
ranges between the purely diffusive case when p = 
(i.e., the Casimir limit) and the purely specular one when 
p = 1. The experimental thermal conductance can then 
be expressed as a function of the effective MFP A e ff. 



K = K C asKff/^Ca 



(3) 



We have normalized the experimental data with respect 
to the Casimir coefficient ficas (according to equation [I] 
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Sample 


L (jLtm) 


Wmin(nm) 


Wmax(nm) 


period (nm) 


MFP (nm) 


Normalized MFP 


b 


P 


straight 


5 


200 


200 





454 


1 


2.27 


0.60 


CI 


4.4 


130 


250 


200 


57 


0.13 


2.52 


-0.32 


C2 


5 


140 


245 


200 


90 


0.2 


2.74 


-0.10 


C3 


5 


160 


260 


200 


178 


0.39 


2.54 


0.22 


C4 


4.4 


230 


250 


50 


70 


0.15 


2.6 


-0.27 



TABLE I: Geometrical parameters and physical data extracted from the thermal conductance measurements of the straight 
and corrugated nanowires (C*). w m in and w max designate the minimum and maximum width of the corrugated nanowires and 
L their length. The thickness of the silicon wires is equal to 200 nm. 



and equation[2| in order to extract the A e ff MFP and de- 
duce the corresponding value of the parameter p for each 
corrugated nanowire. Fig. [3] reports the thermal conduc- 
tance plotted as K/f3cas versus temperature. This nor- 
malization allows a direct comparison between the ther- 
mal transports of each nanowire. From this plot, it is 
already patent that the presence of corrugation strongly 
reduces the thermal transport. 

To go further in the data analysis, a fitting of the data 
K/f3cas versus temperature has been done. From equa- 
tions [I] and [2j the normalized conductance can be ex- 
pressed as 

log(K/0ca,) = a + blog(T), (4) 

where 

& I (>C«s) = \^T b . (5) 

From these two equations, the power law factor b (close to 
3 for 3D materials), the value of the MFP A e jj, and the 
parameter p have been extracted; the data are collected 
in the Table [0 

As it can be seen in this table, the MFPs of the corru- 
gated nanowires are strongly reduced as compared to the 
one of the straight nanowire. The reduction by a factor 
of almost ten is drastic, and can only be attributed to 
the presence of the corrugation. Moreover, for samples 
CI, C2 and C4, the effective MFP is patently smaller 
than the smallest section of the nanowires, demonstrat- 
ing a phonon flux weaker than the one expected in the 
Casimir limit [IB] . Fig. [4] illustrates the different MFPs 
(^straight, A corrug and A Cas ) at the scale of the corru- 
gated nanowires. As demonstrated in reference [9] [10] 
large MFPs in the straight geometry are a signature of 
the presence of ballistic phonons [IT] : those may be par- 
ticularly affected by the corrugation. Finally, the sig- 
nificant reduction in the MFP is also illustrated by the 
p value, which is either small or negative. As proposed 
by Moore et al. in [33], this reduction may be the sig- 
nature of the presence of strong phonon trapping in the 
wire due to multiple phonon scattering or even backscat- 
tering. However, the negative parameter p obtained for 
samples CI, C2 and C4 does not seem to be physically 
sound. 




logT 



FIG. 3: (colour online) Log- log plot of the normalized thermal 
conductance of the straight and four corrugated nanowires 
versus the temperature. The normalization coefficient f3c as 
corresponds to the Casimir limit value of the thermal conduc- 
tance at 1 K. 




FIG. 4: (colour online) The Casimir MFP (A C as = 220nm), 
the effective MFP in the straight nanowire (A str aight = 
454nm) and the corrugated nanowire CI (A corrug = 57nm) 
are compared to the actual geometry of the corrugated 
nanowire. 



In the Ziman model, the paramater p is the probability 
for one phonon reflection to be diffusive or specular. It 
does not take into account the direction of the reflected 
phonons. In our case, as in some recent experiments and 
theorical studies [3 [16 , 20 , 33 , multiple reflections and 
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especially backscattering may have a significant impact 
on the phonon transport. 

If, as suspected, the reduction in the MFP is due to 
the presence of multiple phonon scattering and backscat- 
tering, this should turn up in Monte Carlo simulations. 
Indeed, multiple scattering effects do not involve the wave 
aspect of the phonon but only its corpuscular behaviour, 
something which is particularly well handled by a Monte 
Carlo analysis, which has been carried out here according 
to ray tracing principles [49] . 

In the grey approximation, when the frequency depen- 
dence is discarded, phonons are emitted from one side of 
the wire with a random initial direction and a random 
initial position. The phonons are then tracked until leav- 
ing the wire through one of its sides. The ratio between 
the number of phonons reaching the other side N t hrough 
to the initial number of phonons TV defines the transmis- 
sion T = Nthrough/ N . The only scattering mechanism is 
generated by the surfaces and was modeled by different 
specularity parameters ranging from p = to p = 1 as 
reported in Fig. [5] the value of the parameter p derived 
from the thermal measurement of the straight nanowire 
is p — 0.6. The transmission coefficients can then be 
related to the MFP through the relation jS3 [51] 
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In this equation, A stands for any MFP and L is the 
length of the wire. Following equation [6j a direct cor- 
respondence is obtained between the transmission coef- 
ficient calculated from the Monte Carlo simulations and 
the MFP. The simulations reveal a significant decrease 
in the phonon transmission for the corrugated nanowires, 
corresponding to the backscattering of phonons towards 
the entrance cross section. The comparison between 
the calculated transmissions and experimental MFPs are 
presented in Fig. [5j The simulations provide the cor- 
rect trend when varying the roughness of the nanowires 
through the parameter p: for instance, p is close to 1 
for samples CI, C2 and C4 and a parameter p close to 
has to be chosen for sample C3. As a general conclusion, 
the MC simulations clearly support the experimental fact 
that in the presence of corrugation, the MFP is indeed 
strongly reduced. This reduction may be attributed to 
the occurence of phonon multiple scattering, which is in- 
herent to the ray tracing method. 

Sample C3 presents a smaller mean free path reduc- 
tion: the MFP and the parameter p are smaller than 
that of the straight nanowire but not by as much as in 
the other samples (CI, C2 and C4). Less multiple scat- 
tering can be evidenced in this nanowire. This may be 
due to the nanowire surfaces, which can be rougher than 
in their counterparts, restoring diffusive transport closer 
to the Casimir limit, as confirmed by the MC simulations. 
The behavior of the sample C4 cannot be described by 
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FIG. 5: (colour online) The red full circles are the experimen- 
tal MFP normalized by the MFP in the straight nanowire. 
The Monte Carlo simulations have been performed with dif- 
ferent parameters p (see the figure legend) . The transmissions 
obtained by the MC simulations are normalized to the trans- 
mission in the straight wire. 



the simulations; indeed, the Monte Carlo predictions are 
not able to capture the range of the experimental data. 
However, other mechanisms may be responsible for such 
an MFP reduction. Two scenarios of phonon scattering 
have recently been proposed which could explain a ther- 
mal transport far below the Casimir limit: coherent ef- 
fects, involving correlated multiple scatterings of phonons 
[52] . or strong phonon scattering induced by highly per- 
turbated surfaces (core defects and surface ripples) [53] : 
validating such scenarios would require further experi- 
mental investigations. 

To conclude, the major result of this research re- 
sides in the significant reduction in the phonon MFP in 
monocrystalline silicon nanowires related to the presence 
of multiple reflections induced by corrugated surfaces. As 
a consequence, heat transport has also strongly decreased 
below the Casimir limit. The hypothesis of multiple re- 
flections for phonons is largely confirmed by the ray trac- 
ing analysis showing a strong diminution in the phonon 
transmission. In the best case, the effective MFP has 
been reduced by more than a factor of nine from that of 
the straight nanowire MFP. This is a clear signature that 
phonons are being backscattered due to the presence of 
the corrugated surfaces. This drastic MFP reduction is 
the demonstration that a monocrystalline material can 
be turned into a phonon glass-like material merely by in- 
troducing corrugated surfaces. Corrugations even allow 
breaking the Casimir limit of thermal phonon transport 
in low dimensional systems. The thermal properties of 
these metamaterials is of particular interest in the quest 
for decreasing phonon thermal transport without affect- 
ing the electronic properties, which significantly improves 
the performance of thermoelectrics. 
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